Anharmonicity and self-energy effects of the E 2g phonon in MgB 2 
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We present a Raman scattering study of the E2 3 phonon anharmonicity and of superconductivity induced 
self-energy effects in MgB2 single crystals. We show that anharmonic two phonon decay is mainly responsible 
for the unusually large linewidth of the E^g mode. We observe ~ 2.5% hardening of the i?2 9 phonon frequency 
upon cooling into the superconducting state and estimate the electron-phonon coupling strength associated with 
this renormalization. 

PACS numbers: 74.70.Ad, 74.25.Ha, 74.25.Gz, 78.30.Er 
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High-T c superconductivity in MgB2 is known to be pro- 
moted mainly due to the boron layers 1 , thus the high fre- 
quency lattice vibrations of light boron atoms beneficially in- 
crease the electron-phonon coupling. The E 2g Raman active 
in-plane boron vibrational mode contributes significantly to 
superconductivity; this fact is reflected by the Eliashberg func- 
tion a 2 F(oj) peaking in the same frequency range where a 
high phononic density of states is accounted for by Van Hove 
singularities of the E 2g branch in the T and A points of the 
Brillouin zone (BZ) 2 ' 3 . The reason the E 2g mode plays a 
prominent role in the superconducting (SC) mechanism is that 
the mode strongly couples to the cr-type states of the boron 
plane as can be seen from the basic geometry of the electronic 
configuration^. 

Raman spectra exhibit unusually broad linewidth of the E 2g 
mode 5 ' 6 ' 7 which has been the subject of numerous specula- 
tions. While high impurity scattering in earlier low quality 
samples has been suggested as one of the possible reasons, 
this mechanism can be readily excluded with recent high qual- 
ity single crystals. The two remaining contributions to the 
E 2g phonon rapid decay are (i) strong electron-phonon cou- 
pling and (ii) multiphononic decay (subsequently referred to 
as anharmonicity). The relative importance of the electron- 
phonon coupling and anharmonicity in this matter is under 
intense debate. On one hand a density functional theory cal- 
culation asserts that the anharmonic contribution to the E 2g 
phonon linewidth is negligible (~ 1 me V) 8 . On the other hand 
analysis of the phonon self-energy in the long wavelength 
limit shows that the er-band contribution to the phonon decay 
is vanishing 9 . Thus, even when contributions of the spectral 
weight of a 2 F(cj)| w<W£ , 9 to the damping of the E 2g phonon 

are accounted for, 10 the experimentally observed linewidth of 
25-35 meV at low temperatures 5 ' 6 cannot be explained with 
electron-phonon coupling alone whose part in the E 2g mode 
linewidth at low temperatures amounts to about 6 meV even 
in such an elaborate scenario as that in Ref. 10 . 

Raman scattering experiments have shown that the fre- 
quency of the E 2g mode in single crystals at room tempera- 
ture is around 79 meV 5,6 whereas theoretical calculations sys- 
tematically underestimate this value by about lOmeV- 3 ^. It 
has been suggested that if the E 2g band around the T-point 
is anharmonic then the E 2g mode frequency is increased by 
the missing amount to match the experimentally observed 



valu o 13 ' 14 . In addition, the experimentally observed T c and 
the reduced isotope effect can only be reconciled within 
anisotropic strong coupling theory if the E 2g mode anhar- 
monicity is explicitly included^!!. 

Since the E 2g mode is responsible for a significant part of 
the SC pairing a renormalization of the mode frequency upon 
entering into the SC phase is expected. Moreover, it was pre- 
dicted that the mode renormalization along the T — A direction 
of the BZ is particularly strong 13 . The temperature depen- 
dence of phonon dispersions has been studied by inelastic x- 
ray scattering, however, the expected self-energy effects have 
not been demonstrated^. In earlier Raman work on single 
crystalline MgB2 a change of line width between room tem- 
perature and T c has been observed 6 but not discussed. In the 
experimental context it is important to notice that single crys- 
tal studies should preferably be considered as a valid gauge of 
MgB2 bulk properties because results from powdered, poly- 
crystalline, ceramic, film, wire etc. samples show significant 
variations, especially in relation to phonon excitations (com- 
pare Refsi2«i£). 

In this work we provide spectroscopic information on (i) 
the E 2g phonon excitation resonance, (ii) the reasons for the 
mode's large linewidth, and (iii) the self-energy effects upon 
entering into the SC state. In particular, (i) by means of the 
resonant Raman excitation profile (RREP) we show that the 
light couples to the E 2g phonon via a — ir inter-band transi- 
tion; (ii) we study the E 2g phonon shape as a function of tem- 
perature and fit the linewidth temperature dependence with an 
anharmonic decay model; (iii) we observe the E 2g phonon 
renormalization (mode hardening ~ 2.5%) upon entering into 
the SC phase. 

Raman scattering from the ab surface of MgB2 single crys- 
tals grown as described in Ref. 17 was performed in backseat - 
tering geometry using about 1 mW of incident power focused 
to a 100 x 200 fim spot. We have collected data from two sam- 
ples that are denoted as crystal A and crystal B when a dis- 
tinction is necessary. The magnetic field data were acquired 
with a continuous flow cryostat inserted into the horizontal 
bore of a SC magnet. The sample temperatures quoted have 
been corrected for laser heating. We used the excitation lines 
of a Kr + laser and a triple-grating spectrometer for analysis 
of the scattered light. The data were corrected for the spec- 
tral response of the spectrometer and CCD detector and for 



2 




Raman shift [meV] Excitation [eV] 

FIG. 1: (Color online) Raman response function from crystal A at 8 K for VH (a) and RL (b-d) polarizations corresponding to E2 g channel 
for four excitation energies between 1.65 and 3.05 eV. The spectra exhibit a broad electronic continuum, a superconducting coherence peak 
2 A CT , the E2 g boron stretching phonon mode (/) and two-phonon scattering bands (// - ///). The comparison of a6-plane optical conductivity 
(Ref— ) and the E2 g mode intensity is shown in panel (e). The size of the symbols is proportional to the ratios of mode to electronic background 
intensity. 



optical properties of the material at different wavelengths as 
described in Ref3. 

The point group associated with MgB2 is Dq^; subse- 
quently, polarized Raman scattering from the ab surface cou- 
ples to two symmetry channels: E 2g and A\ g . The corre- 
sponding symmetry channels can be resolved if the appropri- 
ate light polarization is chosen for the in- and out-coming pho- 
tons; in this work only the E 2g symmetry channel will be dis- 
cussed. We denote by (ei„e ouf ) a configuration in which the 
incoming/outgoing photons are polarized along the ei n /e ou t 
directions. In the VH configuration the light is linearly po- 
larized with vertical (V) or horizontal (H) directions cho- 
sen perpendicular or parallel to the crystallographic a-axis. 
The "right-left" (RL) notation refers to circular polarization: 
ej n = (H — iV)/y/2, with e out — e*„ for the RL geometry. 
Both the RL and the VH polarization combinations select the 
E 2g representation. 

In Fig.[T]we show the Raman response as a function of exci- 
tation energy sampling laser lines Xl = 752.5,647.1,482.5, 
and 406.7 nm. We chose steps in excitation wavelengths to 
resolve the broad band at around 2.6 eV in the MgB2 optical 
conductivity^. 

Typical Raman features of E2 g symmetry from a MgB2 sin- 
gle crystal are shown in Fig.Q]c where phonon modes are des- 
ignated with arrows. We see that the E 2g mode at ~ 80meV 
(/.) dominates the Raman response. In addition two-phonon 
scattering peaks 77. and ///. corresponding to flat portions of 
the phonon dispersion curves at half the frequency of the re- 
spective peaks are also observed^. Finally the spectra ex- 
hibit a strong electronic continuum with a SC pair breaking 
peak at 13.6meV corresponding to twice the value of the SC 
gap in the a bands 6 . The evolution of the pair breaking peak 
with cooling below T c is shown in insets of Figs.|2]fl and[3]fl. 

The RREP relates Raman scattering intensity to the inter- 
band transitions seen in the optical conductivity <j(ui). It is 
particularly helpful in the visible range because the interband 
contribution to a a i(u>) contains a pronounced band around 
2.6 eV'^i. In Fig.Q]a-(i we show that the intensity of Ra- 
man response is in resonance with the 2.6 eV band (Fig. []]<?). 



We note that the ratio of the E 2g mode intensity to the strength 
of the electronic background is also largest around the 2.6 eV 
excitation. The Raman intensity of the E 2g mode is plotted 
along the right y-axis in Fig.Q]e as a function of excitation 
while the mode intensity to background ratio is represented 
by the size of the data symbols in the plot. The 2.6 eV peak in 
optical conductivity is associated with the n — > a electronic 
transitions in the vicinity of the T point and a — » 7r transitions 
in the vicinity of the M point of the BZii^. The RREP indi- 
cates that coupling to the E 2g mode is realized via the it <-> a 
interband transition. The strongly resonant behavior of the 
E 2g mode places experimental restrictions on how this mode 
should be investigated. 

In Fig.|2]a Raman spectra for the SC state (zero field) and 
normal state (H > H c2 ) are juxtaposed, demonstrating the 
effect of magnetic field applied parallel to the c-axis at 8 K. 
The E 2g phonon mode hardens below the SC phase transition 
along with a slight renormalization of its linewidth. At the 
same time a SC pair breaking peak at 13.6 meV appears indi- 
cating the presence of the SC gap. In Fig.[2]£> we demonstrate 
how the linewidth of the E 2g mode monotonically narrows 
from room temperature down to T c . To quantify the depen- 
dence of the E 2g phonon frequency and linewidth as a func- 
tion of temperature and magnetic field we fitted the data with a 
sum of phononic oscillators and a SC pair breaking singularity 
on an electronic background as demonstrated in Fig.[2j3. 

In Fig. [3] we evaluate the E 2g phonon frequency u>o{T) and 
the damping constant T(T) for both crystals A and B where 
we distinguish between the respective values for the normal 
and SC states. The solid line in Fig. 3 b is a fit of the damping 
constant T(T) in the normal state to a model of anharmonic 
two and three phonon decay at one-half and one-third frequen- 
cies: 

r(T) = r + r 3 [i + 2?i(r!(T)/2)] + 

r 4 [l + 3n(fi(T)/3) + 3n 2 (fi(T)/3)]. (1) 

Here f2(T) = hcuJh/kBT, with the harmonic frequency 
uJh = 67me V 1 ' 8 ' 23 , n(x) is the Bose-Einstein distribution 
function, Tq is the internal temperature independent linewidth 
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FIG. 2: (Color online) Raman response function in the i?2 9 chan- 
nel for RL polarization and 2.6 eV excitation. Panel (a): Raman 
response from crystal B at 8 K. Spectra at zero field and 7 T are com- 
pared. Spectra are decomposed into two Raman oscillators, elec- 
tronic background (linear slope), and coherence peak singularity for 
the SC state, ~ [u ^/{lo 2 - (2A) 2 )]" 1 . The arrows at the top in- 
dicate the shift of the E^ a mode in the transition from the SC state 
(blue squares) into the normal state (red circles). Inset: evolution of 
low-frequency Raman response with temperature across the SC tran- 
sition. Panel (b): Data from crystal A. Temperature dependence of 
the Ei g phonon mode above T c . Fit parameters to these spectra are 
displayed in Fig.[3] 
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FIG. 3: (Color online) Temperature dependences for the Ei g phonon 
frequency, ujq(T), (a) and damping constants, T(T), (b) for two 
MgB2 single crystals A (black crosses) and B (red circles) for zero 
field cooling (solid symbols) and field cooling at 8 T> H C 2 (empty 
symbols). Solid lines are fits of damping constants T(T) in the nor- 
mal state with eq. (TJ. The parameters from the anharmonic decay 
fits for both crystals are summarized in TableU Inset shows tempera- 
ture dependence of the 2A CT coherence peak which follows from the 
evaluation of spectra shown in the inset to Fig.|2]a. The green bar at 
the y-axis indicates the E^g phonon frequency renormalization upon 
the SC phase transition. 



of the phonon, and 4 are broadening coefficients due to 
the cubic and quartic anharmonicity. The results of the fit to 
this anharmonic decay model are collected in Table Q] For 
both crystals the broadening coefficients T3 + T4 ^ Tq and 
therefore the anharmonic decay is primarily responsible for 
the large damping constant of the Ei g phonon. We identify 
the reason for this rapid phononic decay in the phononic den- 
sity of states (PDOS) peaking at 33 meV, half of the harmonic 
Ei g phonon frequency uj h (see Fig-s,l a in Ref. 3 and Refi2£), 
which corresponds to the Van-Hove singularity of the lower 
acoustic branch (almost dispersionless along the T — K — M 
direction, see Fig. 3 in Ref. 24 ). In this context the narrowing of 
the Eig mode with Al substitution observed in Refs. 5 ' 24 can be 
readily explained with the E-2 g phonon branch moving to ener- 
gies above lOOmeV with increased Al concentration whereas 
the acoustic modes that provide the decay channels stay close 
to their original energies with high PDOS in the energy range 
of 25 - 40meV 24 . In short, the fast decay of the Ei g mode 
is due to the unique combination of its harmonic frequency 



in the T point (67 meV) corresponding to high PDOS at half 
of this frequency. Thus two independent puzzle pieces; (i) 
the cubic anharmonicity contribution being much larger than 
the residual linewidth and (ii) the phonon density peaking just 
at the right energy to provide suitable decay channels, fit to- 
gether, giving evidence that the Ei g phonon anharmonicity is 
responsible for the unusually broad linewidth of this mode. 
In contrast to recent theoretical suggestions 8 current experi- 
ment shows that the Ei g phonon anharmonicity in MgB2 must 
not be neglected. The residual linewidth To that we obtain 
from the fit to the anharmonic decay model, while small, is 
not in contradiction to theoretical estimates 9,10 of the electron- 
phonon decay contribution to the Ei g phonon linewidth. 

It is worth noting that individual Tj parameters differ for 
the two single crystals despite the fact that both samples were 
grown in the same batch. The E% g mode for crystal A is by 
about lOmeV broader than for crystal B. With somewhat 
higher than T® and substantially higher than Tf (see Ta- 
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TABLE I: Comparison of T c and the E2 a oscillator parameters for 
crystals A and B. 



crystal 


T c a 






1 


n 


r fc 


k(%) 


A 


38.2 


79.3 


81.7 


4.0±1.5 


31.4±1.2 


small 


2.8±0.5 


B 


38.5 


78.1 


80.5 


small 


22.9±0.7 


2.8±0.4 


2.3±0.3 



"in K, from SQUID measurements at 2 G; 'in meV 



ble|I|i the E^g mode in crystal A is more anharmonic than in 
crystal B. Accordingly the crystal A mode is pushed to about 
1 .2 me V higher frequency at low temperatures. We note a cor- 
relation between the larger anharmonicity and slightly lower 
T c in the case of crystal A. 

To describe the superconductivity induced self-energy ef- 
fect we refer to Fig.[3]a. Upon cooling in zero field wrj(T) 
exhibits non-mono tonic but smooth behavior down to T c . 
Then at T c it displays abrupt hardening with ojq(T) scal- 
ing to the functional form of the SC gap magnitude 2A CT (T) 
(Fig.0insei). For in-field cooling the E^g phonon frequency 
uiq (T) remains unrenormalized. The differences between the 
phonon frequencies in the normal and SC states at 8 K are 
2.2 ±0.4 and 1.8 ± 0.2 me V for crystals A and B respectively. 
To quantify the relative hardening of the Ei g mode we ob- 



tain the superconductivity induced renormalization constant 
k = (uj$ c /u>o) - 1 w 2.5% (see Table|IJi which is much 
smaller than the theoretically predicted k ~ 12% 13 . 

We estimate the electron-phonon coupling constant 
around the BZ center using approximations adopted in 
Refs^l: A = -KTle(^), where u = n + 
2zcosh~ 1 (o; JV /2A .), and obtain Ag ~ 0.3, which is 
smaller than the theoretically predicted values (compare A up 
to unity in Refs . 1 ' 13 ' 28 ) . 

In summary, we have measured the polarization resolved 
Raman response as a function of excitation energy, tempera- 
ture, and field for MgB2 single crystals. From the temperature 
dependence of the E 2g boron stretching phonon we conclude 
that anharmonic decay is primarily responsible for the anoma- 
lously large damping constant of this mode. For this phonon 
we observe a SC induced self-energy effect and estimate the 
electron-phonon coupling constant. 
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